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Abstract 

Recent work has shown that neutrino oscillations in matter can be greatly enhanced by flips between 
' mass eigenstates if the medium is fluctuating with a period equal to the neutrino oscillation length. Here 

we investigate the effect of the primordial fluctuations on the neutrino oscillations in the early universe. 
We calculate the oscillation probability in the case of a general power law fluctuation spectrum and for a 
' more realistic spectrum predicted by inflation. We also include the effect of the ampliflcation of fluctuations 

resulting from the QCD phase transition. We flnd that there is a region of parameter space where this 
^ • mechanism would be the dominant mechanism for producing sterile neutrinos. However this conclusion 

U ' does not take account of the damping of fluctuations on the neutrino oscillation scale when the neutrinos 

. decouple from the plasma. We flnd that this reduces the probability of flips between the mass eigenstates 

to an unobservable level. 

^ ■ 1 Introduction 

■ For many years there has been tremendous success in determining the masses and mixing of the neutrino sector, 
CO I including experiments on solar [Tj, reactor [HE], atmospheric [4] and accelerator [5l[6] neutrinos. At the Scr level 

the masses and mixing angles determining the solar neutrino oscillations are Am^i = (7-2 — 9.2) x 10~^ eV^ 
\ and sin^ 9x2 = 0.25 — 0.39 [7], while those determining the atmospheric oscillations are Ato|]^ = 2.0 — 3.2 x 10"'^ 

■ and 0.34 < sin^923 < 0.68 with sin^ 2023 < 0.9 [Tj- Currently the sign of Atojs is undetermined and the current 
3(7 upper bound on the allowed values for the final mixing ang le is sin^ 6*13 < 0.044 [TJ. The LSND experiment 
showed oscillations on a mass scale of Am^ ^ ley^[8j if this observation is interpreted as standard neutrino 
oscillations then it can only be explained as oscillations into a forth neutrino state. Precision measurements 
at LEP have shown that there are only three neutrinos that interact via the standard weak interaction exist 

' with mass less than Mz/2 and therefore these addition states must be sterile jF. The global data including 
the LSND (but not including MiniBoone) comprehensively disfavored the inclusion of just one sterile neutrino 
|10] , the inclusion of 2 sterile neutrinos mildly mixed with the 3 active neutrinos has a more acceptable fit to all 
data including LSND [11 although the value of the associated LSND mixing angle is still problematic [12]. The 
inclusion of the MiniBoone experiment did not observe an oscillation [13] but does not ruble out this scheme 
[T4] . Also keV sterile neutrinos which naturally arise in many models [15] can also explain pulsar kicks [16] and 
are warm dark matter candidates . 17] . The latter has recently been received interest following the observations 
of the central cores of low mass galaxies [18] and from the low number of satellites observed in Milky- Way sizes 
galaxies \19^. 

If sterile neutrinos do mix with the standard model neutrinos then they would be produced via oscillations in 
the early universe [20]. If the sterile neutrinos were produced before the active neutrinos decouple (T > MeV) 
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active neutrinos would oscillate into sterile neutrinos, then the active neutrinos would be repopulated, resulting 
in an increase in the number of relativistic degrees of freedom. The expansion rate of the universe would 
increase, which would lead to a higher freeze-out temperature and therefore a higher neutron to proton ratio, 
leading to a higher abundance of helium and other elements [21]. This increased expansion rate can also be 
probed by observations of the CMB [22l [23] . These sterile neutrinos would also contribute to the warm dark 
matter content of the universe which would suppress the formation of large scale structure [23l [24] . 

Neutrino oscillations in the early unverse have been extensively investigated and bounds have been placed 
on the mixing angles and difference in mass squared ^U\. However, these studies have neglected the fact that 
there are temperature fluctuations in the early universe. It has been shown that matter fluctuations can cause 
level crossing between mass eigenstates [3S], resulting in an amplification of these oscillations. In the case of 
active-sterile mixing this could lead to an amplification in the production of sterile neutrinos. If the spectrum 
of the fluctuations was known and the fluctuations were large enough then new bounds could be placed on 
the masses and mixing of the neutrino sector. Alternatively if the masses and mixing was known then then 
constraints could be placed on the spectrum of the fluctuations and on the models which predict them. 

This paper is organized as follows. We begin by calculating the level crossing probability of a two neutrino 
system in an expanding universe with a general spectrum of temperature fluctuations. In sections 3 and 4 
we calculate the level crossing probability for the case of a power law spectrum and a spectrum predicted by 
inflation respectively. In section 5 we investigate the effect of neutrino decoupling and in section 6 the effect of 
the QCD phase transition. Finally we conclude in section 7. 



2 Neutrino evolution 

To consider the Neutrino oscillations in a fluctuating media we follow the analysis developed in PF. The 
evolution of a neutrino system is determined by the Schrodinger equation, i-§fi^ = H*^it)va, where t is time, 
here H'^ and v are the Hamiltonian and neutrino wavefunction in flavour basis. This can be rotated such that 
the hamiltonian is instantaneously diagonalisecH, i^t'm = H^i/'^, where ^'"j = Uu are the instantaneous mass 
eigenstates. For a two neutrino system 

U{t)^( ^^"^ V (1) 



Tj u\ I iddm/dx I , , 



where A„j — y (Am^ cos 20 - Af + (Am2 sin^ 2ef/4:E, tan26l,„ = Sm^ sm2e/{Am^ cos 26 — A), 9 is the vac- 
uum mixing angle, Am^ is the difference in mass squared, E is the energy of the neutrino and A is the difference in 
matter potentials of the two neutrino system. Before the neutrinos have decoupled, for temperatures, T > MeV 
the matter potential takes the form 

A ~ ~Aa,GlE^Tya (3) 

~ -1.86 X 10-20 AC, ( ( .JLJ]\ieV\ (4) 
\MeV J \MeV J ' ^ ^ 



A" 



where a is the fine structure constant and Cy = Un (a — sin^ 6w) sin^ ^?vv/45, a = 3 for and Pe, cl = 1 for 
Dfj^, Vt and v-T and C„ = for sterile neutrinos. We have assumed that the neutrino, anti- neutrino asymmetry 
is of the same order as the baryon asymmetry and therefore the contribution to the matter potential can be 
neglected. There is a resonance when the neutrino system is maximally mixed, when A = Am? cos 26. As 
the matter potential is negative for both neutrinos and anti-neutrinos there will only be a resonance if the 
difference in mass squared is negative. It is convenient to consider the hamiltonian as the sum of the averaged 



^Here we use the notation for a hamiltonian in a medium which is not fluctuating 
^where H^U is diagonal 
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Hamiltonian, and the fluctuated Hamiltonian, 5H, where the true hamiltonian, H = H'^ + SH and in the 
absence of fluctuations 6H = 0. The Schrodinger equation is rotated to diagonahze the averaged Hamiltonian 
where the fluctuations are teated as a perturbation. Now i^i^m — Hm{t)v^i, Hm{t) = H^-^ + SHm where is 
given by Eq. ([2]) and 



SV ( -cos20„ 



sin 26r, 
cos 29. „ 



(5) 



where V = A/2E and 5V = V ~ V[), where Vb = Ao/2i? and Aq is the matter potential in the absence of 
fluctuations. If the off diagonal elements are non-zero transitions between the mass eigenstates can occur. If 
Am ^ \d9m/dt\ the neutrino propagation is non-adiabatic with or without fluctuations. If \d9m/dt\ <^ A™ 
the evolution in the absence of fluctuations is adiabatic and the level crossing between the mass eigenstates is 
determined by the off diagonal elements in Eq which are 5V sa\29 /2. In the perturbative limit where the 
fluctuations are small the level crossing probability is 



P: 



at exp 



dt'2l^r, 



(6) 



Numerical simulations have shown that for large fluctuations the system becomes depolarized, i.e. with the 
probability of detecting each flavour of neutrino being 1/2; a rough criteria for depolarization is the region 
P > 1/2. The neutrinos that we consider are in thermal equilibrium with the plasma and therefore have an 
energy spectrum which is approximately Fermi-Dirac and is charectorised by the temperature of the plasma. 
For this spectrum the average energy of the neutrinos (E) — yT, where y ~ 3.1514. For the rest of this 
paper we consider the oscillations of a neutrino with an energy equal to the average energy of the ensemble of 
neutrinos, i.e. E = yT. To calculate the fluctuated hamiltonian for a neutrino with this energy we use Eq.([3]) 
V = A/2E (X ET'^, in the case of smaU fluctuations 6V = W5T/T + V6E/E. If the neutrino is in thermal 
equilibrium with the fluctuations SE/E = ST/T and therefore SV = 5VST/T. If the neutrino is not in thermal 
equilibrium with the fluctuations then SE/E ~ 0, i.e. the energy of the neutrino remains approximately constant 
across the scale of the fluctuation. As we will see in the section [5] the neutrinos decouple from the fluctuations 
at the scale which amplifies the neutrino oscillations before the time at which the neutrinos oscillate. Therefore 
it is a good approximation to use assume that the neutrinos are not in thermal equilibrium with the fluctuations 
and SV = 4yS. Expanding the fluctuations into their Fourier components this now becomes 



P = 



(2^ 



Vik) 



dt 2V sin 29„i exp 



dt' i 2Am 

a 



(7) 



where k is the comoving wavevector, 'P{k) = |i5T(fc)p, Srik) is the Fourier transform of the temperature 
fluctuation which we have assumed to be constant in time, a is the scale factor and z is defined to be the 
direction of neutrino propagation. In this derivation we have approximated the neutrino to be massless, therefore 
ds^ ~ dt^ — a{t)'^dx'^ = 0, and |Ax| = J dt/a. Note that the integral in Eq. ([7]) is oscillatory except when 
2Am = kz/a, i.e. the fluctuation length is equal to the neutrino oscillation length. To solve the time integral 
we use the stationary phase approximation. Here we define new variables 



Q{a,kz) 
q{a) 

F{a) 
A 



q{a) - kz/h, 

— A„i = ^ (a + Aa j , 



2Va sin 29 „ 



2yToh' 



2/iAsin26', 



2AC^GlyTl 
a Am? cos 29 ' 



(8) 
(9) 

(10) 

(11) 

(12) 
(13) 
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where for temperatures of interest the universe is in the radiation dominated era, where T = To/a, y = E/T, the 
hubble parameter H = a/a = h/a^ and Tq, y and h are constant. We define |/p such that P = J (^J^i 'Pik) , 
applying the stationary phase approximation 

2 

daF{a)e'J''^'^'''''''^'^''' 



2TrF{a,i{k,))^ 



27r 



+ (a. 



as2, 



(14) 
(15) 



(asi as2) 



\Q'{a,iik,),h)\ 

27rf(a,i(fc,))F(a,2(fc^)) 
\Q'{a,,{k,),k,)Q'{as2ik,),k,)\ 



asl(fez) 



(5(a, kz)da. 



(16) 



aai(fe^) 

The integral in Eq. is oscillatory except where Q{k,a) = 0, this defines o^i where i=l, 2, such that 

Q(asi,kz{asi)) = 0, note there are two real solutions to this equation. Once |/p is integrated with respect to 
kz the final term of Eq. pop is averaged to zero and therefore it is no longer considered. Integrating |/p with 
respect to kz is done so by a change of variables, dkz = ■^^da^i = hQ' {asi)dssi- Now the limits of integration 
are from kz = [—00,00], this corresponds to a^i = [0, oq] and a^i = [ao,oo], where oq. Therefore the integral 
can be simplified to 



P= — I dkr 

27r 



da krF{afV{^/kl + kz{aY), 



(17) 



where we have introduced the cylindrical polar co-ordinates (fc^, kz,9). Note for a^i there is a minus sign from 
Q' {asi{kz), kz) = —\Q'{asi{kz),kz)\ and also the limits of integration, where kmin = and k^ax = 0. This 
results in a positive term in Eq (jl7p . Now to calculate the level crossing probability we need the spectrum of 
the fluctuations. 



3 Power Law Spectrum 

The measured spectrum of fluctuations is approximately fiat. This correlates to 7^ = 27r^A:~'^A^(fc), where 
A^(fc) — A^(/co)(^/^o)"°~^i and the spectral index — 1 is small. From the combination of the WMAP3 [55] 
and SDSS [57] data, Ug — 0.98 ± 0.02 at the 68% confidence level |5S]. To first approximation we assume that 
the spectral index is constant. Our motivation for this is to analyze the effect of how the parameters effect the 
level crossing probability for the simplest case possible. For this case the level crossing probability is 



A^iko] 



//i^cos26iAi/3^ 













cos 26* 



sin^ 26*- 



(l+u6)3-(n=-l) 



du, 



(18) 



where we have defined u = X^^^^a. The first square bracket is the size of the fluctuations at the neutrino 
oscillation scale with a scaling factor of A^^'^. From Eq. ([7]) we can see that the fluctuations at the neutrino 
oscillation length scale is the important factor in increasing the level crossing probability. The second square 
bracket is the weighted number of oscillation lengths. This can be more easily seen if one expands fiX^^^ — 
/{2yTo/ Am'^X^/'^) and realising that is the scaled Hubble length and (2yTo/Am^A^/^) is the scaled 
oscillation length. As one might naively expect the level crossing probability increases with the number of 
oscillation lengths. There is also a mixing factor sin^ 29 and a numerical factor which are as one would expect. 
Considering the oscillation of electron neutrinos into sterile neutrinos the oscillation probability is given by 



Pes = 1/2 - (1/2 - F) cos 26*, cos 26 



(19) 



where 9i and 9f are the mixing angles at production and detection of the neutrino. For small mixing cos 2^ 
1 — 29"^ and therefore 



{I - 2P){0^ + 0^f) + P. 



(20) 
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For > 1 the oscillation probability is dominated by the level crossing probability caused by the fluctuating 
plasma and is this dominant mechanism for the production of sterile neutrinos. In Figure [1] we plot the contours 
of constant P/0^ as as function of mass squared differences and spectral gradient. The level crossing probability 
increases with Am^ , this is because this corresponds to smaller oscillation length and therefore more oscillation 



periods. Also the level crossing probability increases with 
on the neutrino oscillation length scale are larger. 



this is because the amplitude of the fluctuations 
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Figure 1: The contours of constant as a function of spectral index 71^ — 1 and difference in mass squared 

Ato^. 



4 The primordial spectrum from Inflation 

An initial stage of inflation seems a necessary ingredient ot the Big Bang model [29]. An intriguing feature of 
inflationary models is that they predict a nearly scale invariant spectrum which is consistent with that observed 
by WMAP 26J, the 2dFGRS [30] and SDSS [27]. A simple model capable of generating this spectrum of 
fluctuations for inflation is a single scaler field, 0, with the potential V{(j)), in the slow role approximation 
e = {MpiV /V)'^/2 < 1 and \ri\ = M'^i\V" /V\ < 1, where ' represents d/dcj). For this model the curvature 
fluctuations, A^(fc), which are related to the temperature fluctuations by 'P{k) = 27r^fc^^A|j(/c) are given by 



A?j(fc) 



i/4 



(27r0)2 



(21) 



where <, is the time at which the k mode crossed the horizon, k = aH\t^. This together with the equation 
of motion for the inflaton field in the slow role approximation, iHcj) — —V , and the Hubble parameter, 
ij2 = SttV/SM^i resuhs in 



Approximating H to be constant and using k = aH\t, 

ln(fc/fco) = - 



(22) 



irVdct) 



pi 



(23) 
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For simplicity we consider the potential V = Uncj)"', for this case 



pi 

where Vb = UnipQ- We can see that the spectrum is flat for n = 0,-2, i.e. for Ug = 1, increases with k for 
— 2 < A: < 0, for > 1 and decreases with k for n > and n < — 2, for < 1. Now to obtain analytical results 
we can use observations of the primordial spectrum to determine the three of the four free parameters. At the 
68% confidence level In (A|i,(fco)10i'^) = 3.17 ±0.06 and dlnAydlnk\k„ = ns(fco) - 1 = -0.02 ± 0.02, where 
ko = 0.002Mpc~"'^ [28j, re-expressing Eq ((24|) in terms of these parameters 



A^(fc) = A^(fco)ll+ ''^^^°^ ^ ln(fc/fco)) , (25) 



where we have chosen a„ = (2 + n)/2 to be undeterminecjf|. 

For this spectrum of primordial fluctuations the level crossing probability is 



P = 27rA^/^Ai/2 sin^ 261 / dv du^ (26) 

Jo Jo 

/ 2 \ — 3/2 

J. ^ ^ ^ (l - (in {v + {u^ + u-r) + ln(VA^/Vfco)))"" (27) 

Figure [4] plots the contours of constant P/d^ for n ~ —1 for rig > 1 and n = 4 for rij < 1. The fluctuations 
increase with Am^ and Ug for as for the case of the power law spectrum. However the increase with Ug is 
not as pronounced as there is only a logarithmic increase with Ug. This mechanism becomes dominant in the 
production of sterile neutrinos for keV neutrinos with > 1. 




Figure 2: The contours of constant level crossing probability for a spectrum predicted by inflation with a 
potential V — a„(/)" where n = — 1 for Ug > I and n 4 for < 1. 



^For fluctuations to be present at the neutrino oscillation length scale we require that the number of e-folds, A'^e = J da/a ~ 
60 < ln(feo/fc), for mass squared differences up to keV^ this is true. 



6 



5 Neutrino Damping 



So far we have neglected the damping of the fluctuations in the early universe by neutrino diffusion which 
is equivalent to Silk damping of photon decoupling. This occurs when the neutrinos diffuse from over dense 
regions dragging the matter with them. Neutrinos decouple at T^^'^ ~ lAMeV for i^e and T^^^^^ — 2.2MeV for 
and Vr |31| . this occurs when universe expands faster than the neutrino interaction rate H > T^, where Ti, 
is the neutrino interaction rate. However, the neutrinos decouple from a particular mode of oscillation when 
Cskph > Fi/, where Cs is the sound speed of the oscillation. Following |32j the damping factor in the density 
fluctuations is 

D{k, rf) = exp (-F (fc, 77)) = exp j {kphys/etot)ilvisckdri^ (28) 

where ptot is the total energy density, %isc ~ (4/15)S/9^„Tjy^ is the shear viscosity co-efficient, p^^ is the energy 
density of Va, t^^ = 1/Fi,^ is the typical collision time and T^^ is the interaction rate of Va with the medium, 
Vmax — "im[77, 77dec(fc)], Vdecik) IS the comoviug time at the decoupling of the oscillation mode k and the sum 
is over all active neutrino species. Substituting the neutrino interaction rate Fi. = ^^G'^T^, the energy density 
ea — 7r^(7a(T)r'*/30, where g is the effective number of relativistic helicity degrees of freedorr0 into Eq. ([28|) 

where Tmax = max[T, Tj^'^'^(fc)], Tj^'^'^(fc) decoupling temperature of ly with the oscillation mode k. The decoupling 
temperature of an oscillation mode can be related to the neutrino decoupling temperature by 

rf^(fc)=Tf=(^^)^|^ , (30) 

where T^'^'^ is the ly' decoupling temperature. For all k modes the neutrino has decoupled before the neutrino 
oscillation length is equal to the oscillation length (see appendix A) and therefore T^ax = Tdedk) and the 
damping is maximal. This damping modifies the primordial spectrum to 

A2(fc) ^ A2(fc)exp(-F(fc)), (31) 



T{k) = A(^{u^ + u-''f+v^ , (32) 

where 

A-E4— (t^) c-^^'U.eXilcos2eY'\ (33) 

where Xdec = ^oT^tc^ f-'-dec — I^O'dec '^dec is the scale factor at the time of neutrino decoupling. This effectively 

gives a cutoff in the fluctuation spectrum when + + — A^^^^ — 3.5 x 10~^° (sm^) ■ 

A cx A (m^) the cutoff is larger for smaller mass differences. This is because smaller mass differences 
corresponds to a larger oscillation length, therefore fluctuations with smaller k. From Eq. ((29l) we see that 
damping is largest for (a)larger fc/F ^ t/Losc, where Lose is the oscillation wavelength and {h)k/H. Physically 
(a) is due to the neutrinos 'dragging' the matter to distances larger than the oscillation wavelength and (b) is 
for smaller expansion rate and larger oscillation frequency there will be a larger number of oscillation lengths 
when the damping is large. The effect of this damping on the power spectrum with ris = 1 is shown in Figure [31 
Here we see that the damping severely reduces the level crossing probability for Am^ > 10"''. This damps any 
oscillations into sterile neutrinos in the eV — keV mass range and the oscillations between active neutrinos. 

*where fermionic degrees of freedom are suppressed by a factor of 7/8 
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Figure 3: The effects of damping on the power law spectrum with and without fluctuation caused by the QCD 
phase transition. 



6 QCD phase transition 

At a temperature of about T* ^ IbOMeV there is a QCD phase transition from a quark-gluon plasma to a 
hadron gas, this transition can increase the amplitude of the primordial fluctuations for sub horizon modes 
while leaving the superhorizon modes unaffected [33j . According to the second law of thermodynamics 

P + P = T^ (34) 

where p and p are the energy density and pressure respectively. During the phase transition the energy density 
is discontinuous in temperature at the temperature of the phase transition T*, therefore the pressure must be 
continuous with a discontinuous gradient. During the phase transition the energy density decreases from p_(r*) 
to /9_|_(r,), where we now use the notation — and + for the beginning and end of the transition respectively, the 

1/2 

pressure p remains constant plT^), and therefore the speed of sound Cg = (dp/dp) J = 0. As the sound speed 
is zero there are no restoring forces from pressure gradients and therefore the radiation fluid goes into free fall. 
This can be seen quantitatively from the equations of motion for sub horizon density fluctuations 

<5' - fc^A = (35) 
-ip' + cl5 = Q (36) 

where ip — {p/ {p + p))vpec and Upec is the peculiar velocity. During the phase transition ip' ^ and S' = 
kip = const. The solution to these equations are — ip- and S = + k{r/ — 77_)i/)_, i.e. the energy density 
increases linearly with the time of the phase transition and the wavevector. This results in the ratio of the initial 
fluctuations, Ai, to the final oscillations, Af being Af /Ai ~ k/ki where ki = y/3/{r]+ — ri-). This approximately 
modifies the spectrum of density fluctuations by 

a2, ,^ f A2(77,fc) fork^ki 
^ ^'''''^ \ A2(,7,fc)(fc/fci)2 fork>ki 

P/O"^ is plotted in Figure [3] including the effects of the QCD phase transition. From this we see that the level 
crossing probability is increased but the fluctuations are still severely damped and consequently there is still no 
observable effect. 



7 Conclusion 

If sterile neutrinos mix with the standard model neutrinos they will be produced by oscillations in the early 
universe and the number produced will depend on the mass squared differences and mixing angles. Sterile 
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neutrinos would add to the relativistic energy density of the universe increasing the expansion rate of the 
universe resulting in an increase in the abundance of light elements. They would also contribute to the warm 
dark matter of the universe. Neutrino oscillations propagating in a fluctuating medium can be enhanced due 
to the increase in the level crossing probability between mass eigenstates. This probability is dominated by 
the amplitude of the fluctuations on the neutrino oscillation scale. The temperature fluctuations in the the 
early universe have been probed by the WMAP, 2dFGRS and SDSS data and have shown that this spectrum of 
fluctuations is nearly scale invariant. In this paper we have calculated the level crossing probability for the case 
the spectrum of density fluctuations has a simple power law dependence and for the case this is modified by 
logarithmic running as predicted by a simple model of inflation. For these cases there is a region of parameter 
space for large Am^ and rig were the production of sterile neutrinos would be dominated by this mechanism. 
However this conclusion ignores the damping of the fluctuation which occurs when the neutrinos decouple from 
the plasma, this is equivalent to Silk damping. We have shown that the neutrinos damp the fluctuations on 
the scale of the neutrino oscillation length scale and thus the level crossing probability is vastly reduced and is 
unobservable. A further effect is that the fluctuations on the neutrino oscillation scale are amplified at the QCD 
phase transition. However we find that this enhancement is insufficient to overcome the damping that occurs. 
In conc;Iusion the; production of sterile neutrinos docs not get amplified from level crossing between the mass 
eigenstates due to primordial fiuctuations because when the neutrinos decouple they damp the fluctuations at 
the neutrino oscillation scale. 
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Appendix A 

The decoupling temperature of the k oscillation mode is 



, X 1/4 

^(fc) = r,^- ( ^ j . (38) 



The wavevector of the the neutrinos oscillation is 

k = hn\^'^ cos 29^ V + {v? + u-'^f (39) 

where u = X^^^^a = X^^^^T/Tq where T is the temperature of the universe. The ratio of the decoupling 
temperature to the temperature of the universe is 



= f Cs^ldecXdec COS 26u^\Jv+ {v? + U-'^f \ 



(40) 



where A^ec = ^^dec^ l^dec = Mo^dec ^'^'^ '^dec is the scale factor at neutrino decoupling. The right hand side is 
minimized for v = and u = 0, where at this minimum 



= ((3/22/3)c,MdecAdec COS 2^) ~ 2.5. (41) 

As the minimum of this ratio is greater than unity the temperature of the universe when the wavelength of 
the fluctuation is equal to the neutrino oscillation length is always greater than the temperature at when the 
neutrinos decoupled from this mode. 
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